WUE-ITP-96-006 



hep-ph/9605420 



High-Energy Approximation of One-Loop Feynman 

Integrals 

M. Roth and A. Denner 

Institut fur Theoretische Physik, Universitdt Wiirzburg 
Am Hubland, D-97074 Wiirzburg, Germany 



Abstract: 

We provide high-energy approximations for all one-loop scalar 3- and 4-point functions 
and the corresponding tensor integrals that appear in scattering processes with four ex- 
ternal on-shell particles. Our expressions are valid if all kinematical invariants are much 
larger than the internal and external masses. They contain all leading-order terms of the 
integrals. 



WUE-ITP-96-006 
May 1996 



1 Introduction 

The experiments at LEP and the SLC have provided a large number of high-precision 
experimental data. In order to compare these data with theoretical predictions the inclu- 
sion of one-loop radiative corrections has been inevitable, and even the leading two-loop 
corrections had to be taken into account [ |TJ. Although forthcoming experiments as e.g. at 
LEP2 or a Next Linear Collider (NLC) will not reach such a high experimental accuracy, 
nevertheless the inclusion of one-loop radiative corrections will in general be required in 
the corresponding theoretical predictions. In fact, as the radiative corrections typically 
involve logarithms of the energy, their relevance even grows with increasing energy. 

The radiative corrections to e + e~ —►//(/ ^ t) are relatively simple even off the Z res- 
onance, because only external fermions are involved, and their masses can be neglected at 
sufficiently high energies. Once the masses of the external particles are relevant, as e.g. 
for top-pair production, or external gauge bosons are present, the number of Feynman 
diagrams is considerably increased, and the evaluation as well as the results for the one- 
loop radiative corrections are much more complicated. In addition, owing to the presence 
of different energy scales, numerical instabilities eventually arise in the TeV range or be- 
yond when the standard calculational procedures are used. These instabilities are caused 
by numerical cancellations originating from the recursive reduction of the tensor inte- 
grals and — in reactions involving longitudinal gauge bosons — additionally from unitarity 
cancellations. Evidently, it is desirable to circumvent these problems. 

As the energy scale for future experiments is large compared to the electroweak scale, 
a natural approach consists in neglecting all masses as compared to the centre-of-mass 
energy whenever possible. Such an approximation has been worked out for the process 
e + e~ — > W + W~ [ 0. This calculation has shown that the one-loop corrections indeed 
simplify a lot in a high-energy approximation, and that such an approximation can already 
be useful for energies around 500 GeV. 

With purely algebraic manipulations, which can easily be performed by computer- 
algebra packages such as e.g. FeynCalc [ |3|], all one-loop Feynman amplitudes can be 
expressed in terms of scalar and tensor integrals. The tensor integrals can be algebraically 
reduced to scalar integrals [ [|. Moreover, all scalar integrals with more than four internal 
propagators can be related to scalar 4-point functions [ |5|. Consequently, all one-loop 
Feynman amplitudes can be reduced to scalar 1-, 2-, 3-, and 4-point functions. For these 
functions complete analytical results exist [ [|. 

Thus, the expansions of the scalar integrals for high energies are an essential ingredient 
of a high-energy approximation. The example of W-pair production has shown that the 
evaluation of these approximations requires a substantial amount of work. In order to 
facilitate future calculations in the high-energy limit, it is therefore desirable to calculate 
and tabulate the high-energy expansion of these integrals. 

In the literature, methods have been described to evaluate Feynman integrals in the 
limit of large masses and/or momenta [ [7]. However, these methods in general do not 
apply when the external particles are on their mass shell [ BJ. In this case (IR or mass) 
singularities can arise if all masses are neglected. For UV-singular Feynman integrals, 
including all 2-point integrals, no IR or mass singularities show up such that the corre- 



1 



sponding high-energy approximations are simply obtained by putting all masses to zero 
in the exact results. 

One can, of course, use the known exact results for the 3- and 4-point functions to 
derive the approximations. However, the necessity to keep the masses finite in order 
to extract the correct singularities renders the expansion of the exact results a tedious 
exercise. 

On the other hand, the Mellin-transform technique allows to construct approximations 
in a relatively simple and direct way as long as only one of the kinematical invariants gets 
large, as it is the case for 3-point functions. Moreover, it turns out that this approach can 
also be used for the scalar 4-point-function despite of the fact that two invariants become 
large in the corresponding high-energy limit. However, the Mellin-transform technique 
fails for the tensor 4-point integrals. Therefore, we used the tensor-integral reduction 
algorithm to construct the corresponding high-energy approximations. In this way we 
have calculated high-energy approximations of the basic one-loop integrals necessary for 
scattering processes of two on-shell particles into two on-shell particles, i.e. for the most 
important processes at high energy colliders. 

The paper is organized as follows: In Sect. |2| we describe our evaluation of the high- 
energy approximation for the scalar and tensor 3-point functions and list the correspond- 
ing results. In Sect. |3| the same is done for the 4-point functions. In the appendix we 
summarize our conventions for the Feynman integrals. 



2 High-energy approximation of one-loop 3-point functions 

The scalar and tensor 3-point functions are defined in App. |A.1| . We first describe the 
evaluation of their high-energy expansions and then list the results. 

2.1 Calculation 

In order to calculate the high-energy approximation of the scalar 3-point function, we 
start from the Feynman-parameter representation 



n ( \ AAA 5{l - X - X X - X 2 ) 

C (pi,P2,mo,m 1 ,m 2 ) = - dx dx 1 dx 2 — / r -, (1) 

Jo g(x , xi,x 2 ) - (pi - p 2 ) 2 x 1 X2 

with 



g(x , x l: x 2 ) = mlx + m\x x + m\x2 - p\x a xi - p\x x 2 - ie. (2) 
We search for an asymptotic expansion in the limit 

r = \(Pi ~P2) 2 \ > rnl,ml,ml,pl,pl, (3) 

including all terms of order l/(pi — P2) 2 ■ Looking at the specific results derived for the 
high-energy approximation of the W-pair production [ §], we notice that the scalar 3-point 
function has the asymptotic form: 



C (p 1 ,p2,m ,m 1 ,m2) = J (^ln 2 (r) + c x ln(r) + c ^j +O (-^) . r — oc. M i 
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This asymptotic expansion cannot be extracted by performing a simple Taylor expansion 
or by using the general methods described in Ref . [ [7|] . It is convenient to use the Mellin- 
transform technique, because each pole of the Mellin transform is uniquely related to one 
of the terms in the asymptotic expansion (f|). 

If the function f(r) fulfills certain integrability conditions (see Ref. [0), the Mellin 
transform _ 

M\f,£]= / r«"V(r)dr (5) 



converges absolutely and is holomorphic in a vertical strip a < Re(£) < (3 of the complex 
plane. The corresponding analytical continuation contains poles in the complex half- 
planes Re(£) > (3 and Re(£) < a. The inversion of the Mellin transform 

1 rc+ioo 

f(r) = — r-*M[f,m with a<c<(3 (6) 

Z7T1 J c— ioo 

can be performed by closing the integration contour over the right-hand complex half- 
plane at infinity. According to residue theorem f(r) is obtained from the poles of the 
Mellin transform in the right-hand half-plane Re(£) > (3. 

The inverse Mellin transform of a single pole of order (n + 1) is given by 

— / d£r-*- — — = ^ln», (7) 

i.e. a term that appears in the asymptotic expansion (|j). From this it is obvious that 
we get an asymptotic expansion for r — > oo by closing the integration contour in the 
right-hand complex half-plane. Moreover, the leading terms originate from the poles 
of the Mellin transform A4[f, £] on the right edge of the convergence domain, i.e. with 
R- e (£o) = P- Thus, we can write down the poles of the Mellin transform that determine 
the asymptotic expansion (|4]) of the scalar 3-point function 

M[C , = + ir ^ ? + ^ + 0((1 - 0°) (8) 

with the same constants c$ as in (U) [see (0)]. 

In order to evaluate (multiple) Mellin transforms of one-loop integrals we use the 
generalized Feynman-parameter representation [ [[(J 

where £ = X)£=o&! R- e (&) > 0, «j > 0, and the A4 are not on the negative real axis. 
Choosing «j = for % > and «o = 1 and performing the integration over r we get the 
following formula for a multiple Mellin transform: 

dri---dr B71 ^ n A r ? = (i ) 

with £ = £ — EILi Ci- This formula can directly be used to evaluate the Mellin transform 
of ([]]) upon identifying: 

A = g(x ,x 1 ,x 2 ), At — ( <j '\t)X\X<i and r x = r = | (pi - p 2 ) 2 \ (11) 
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with a = (pi —p 2 ) 2 /r — ±1 and n = 1. In order to ensure the validity of (10), an infinites- 
imal imaginary part of the ie-prescription is given to both A and A\. After integration 
over Xq we find the Mellin transform of the 3-point function, which is holomorphic in the 
vertical strip < Re(£) < 1, 

M[C ,£] = - U 1 , / dx, / dx 2 - — g-j (12) 

f-cr — ie)« Jo Jo h{x ll x 2 ) l ~ i x\x\ 



with h(xi, X2) = g(l — X\ — X2, X\, x 2 ). According to the above discussion, the asymptotic 
expansion is determined from the poles at £ = 1. As T(l — £) includes a single pole, we 
need the expansion of the integral about £ = 1 including the order 0((1 — 

First we assume h(0, 0) = ?71q 7^ 0. Naive expansion of h(x 1 , x 2 )~^ 1 ~® about x\ = 
and x 2 = would yield two infinite series of pole terms contributing to the high- 
energy approximation namely the Taylor expansions of h(xi, O) - ^ 1 "^ and h(0, x 2 )~^~^\ 
Therefore, we use 

1 ' ' 1 — + 0(x 1 x 2 )0(l-£), (13) 



h(x 1 ,x 2 ) 1 -Z hix^oy-t h(0,x 2 ) 1 -t h(0,0) 

such that the terms that are not explicitly written out involve the product of X\ and x 2 , 
and get a compact expansion of the Mellin transform: 

m\c - r^-fl f 1 

M[C °^ ~ ~ { - a - i€) t { m o)w ( i _ 2 - y 

1 dx- In f M - t dxi In fM) + Oil - £) } . (14) 

x v>(o,o)y io x \h{0,0) J J 

From the inversion formula @ or by comparing (f|) with (§) we obtain the high-energy 
approximation of the 3-point function: 

_ /' dl I ln fM) _ /' dl I h (Ml + o (. 1 J . (15) 



a; VMM)/ 70 x VMM)/ J V(pi-^) 4 

For h(0, 0) = = the expansion ( |T3"D breaks down. Therefore, we extend the 
integration domain of the Mellin transform (|12|) to a square (0 < Xj < l,i — 1,2) and 
divide this square into two sectors^: 

x\ > x 2 and x 2 > x\. (16) 

It can easily be seen that the extension has no poles at £ = 1, apart from the one contained 
in the pref actor T(l — £), and the integral yields 

1 dx 1 f 1 dx 2 - 1^^ = ^ + 0(1-0. (17) 

Jl-xi h{XiiX2) K x\x 2 " 



1 This sector decomposition is similar to the one used in the proof of the convergence theorem for 
Feynman amplitudes [ O] . 



By substituting x 2 — > X\X 2 in the first sector and X\ — > x 2 X\ in the second sector we get 
the following result for the Mellin transform: 

M[C ,Z] = - r | e) 5 (1 ~f { t^2— _ r, n 



— cr - le 



o Jo (h(xi,xix 2 )/xi) x\xl 



+ C^I>* 1m^^M -^ + 0(1 -^ (18) 



with 



h(x!,xix 2 )/xi = m\—p\ + (m\ - p\)x 2 + (p\ + p\x 2 )(l + x 2 )xi - ie, 
h(x 2 xi,x 2 )/x 2 = m 2 2 -p 2 2 + (m\ — pQxi + (pjxi + pl)(x 1 + l)x 2 - ie. (19) 

For m\ 7^ p\ and m\ 7^ p\ the two integrals in (|T8| ) can be evaluated similarly to ( |T2| ) in 
the case 7^ 0. 

If m\ = p\ 7^ we perform a further sector decomposition in the first integral of ([18]). 
The case m 2 = p\ 7^ is covered by exchanging the indices 1 and 2. The remaining cases 
are either IR-singular (m = 0, m\ = pf, m\ = p 2 ) or mass-singular (wq = m\ — p\ — 
or mQ = mj = P2 = 0)- Thus, we can construct high-energy approximations for all scalar 
3-point functions that are neither IR- nor mass-singular. The singular integrals are also 
covered by our results if the singularities are regularized with masses. Our general results 
were checked by comparing with the specific expressions in Ref. [ |2|]. 

The high-energy approximation of tensor 3-point integrals can be evaluated in the 
same way: 

Consider first the coefficient functions C\...\. They have the same Feynman-parameter 
representation as the scalar 3-point function ([!]), but with an additional factor x\ in the 
numerator [see (|A20 )] which prevents the Mellin transform from being singular at £1 = 0. 



Consequently, the integrand can be expanded with respect to x 2 , and terms of order 
0(2:2) are not contributing to the high-energy behaviour. As a consequence, the results 
are independent of m\ and p\. 

The functions C\...\ 2 ... 2 have no singularities at x\ = and x 2 = owing to the 
numerator factor x\x 2 . The function T(l — £) contributes the single pole 1/(1 — £). Since 
we can set £ = 1 in the integral, the result is independent of the internal and external 
masses. 

The remaining coefficient functions of the 2- and 3-point functions are UV- divergent. 
They have no IR or mass singularities, because the integrand of the Feynman-parameter 
representation has a logarithmic structure. It can be easily checked with the help of the 
Mellin-transform technique that the high-energy approximation is obtained by neglecting 
the internal and external masses. 

In contrast to the asymptotic expansion of the scalar 3-point function, those of the 
coefficient functions of the tensor 2- and 3-point integrals include no dilogarithms and 
need not be split into several different cases. All results for these tensor integrals were 
derived using the Mellin-transform technique and checked by tensor-integral reduction. 

2.2 Results 

The following results hold for 3-point functions in the limit 

\(Pi -P2) 2 \ > rnl,ml,ml,pl,pl, (20) 
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and for 2-point functions in the limit 

\p 2 \ >m 2 ) ,m 2 . (21) 

2.2.1 Scalar 3-point function 

For the scalar 3-point function we have to distinguish several different cases: 

ml ^ 0: 

n , \ 1 J 1 i 2 (-(Pi -P2) 2 — ie\ 

(Pi-P2) 2 l 2 V m o- lf / 

+ m 0) mi) + Ic(P2> w o,m 2 ) j, (22) 

where 



r 1 1 / m? — m 2 , — p? p? 
IciPi, m o, m i) = — I dx- In ( H ^ : x H 2 — r x 



2 

<; 

x \ tuq — le wig — le 

= Li 2 ( — 7) 9 5 — ; 7—5 5 ; 5 — r ) (23) 

± ymQ — mf + pf ± K{pf, — le, mf — le) j 

with /t(a, 6, c) = Vet 2 + & 2 + c 2 — 2a6 — lac — 2bc; 

ml = 0,m 2 ^ p\,m\ ^ p\: 

nt \ 1 f 1 f -(Pi ~P2) 2 -i^ , / -(gi -P2) 2 - ie \ , 

Co(pi,P2,m ,mi,m 2 ) ~ 7 rj^ In g 5 : — In : — (24) 

[Pi ~ P2) 2 { \ mf — pf — le J \ mi -pi -is 



2 2 1 z l 2 2- 

mf — pf — is J \ mi — pi — ie 



m 2 = Qj m 2 = p 2 ^ Qj m 2 ^ p2. 



, 1 f, (-(P1-P2) - ie \, (-(JP1-P2) -ie\ 

Go{pi,P2,mo,m 1 ,m 2 ) ~ rr< In = = — : — In = — : (25) 

{Pi-P2r{ \ mi -pi -ie \ pf-ie 



A.( -^-ri- : ^ J-^)^|. 

2 V m 2 - P2 - ^ / V m 2 - P2 - ^ / 6 J 

The remaining cases are found by exchanging the indices 1 and 2 or are IR- or mass- 
singular. Specific examples of these general results can be found in App. A of the second 
paper of Ref. [ 

2.2.2 Tensor 2- and 3-point functions 



The coefficient functions of the tensor 2- and 3-point integrals, defined in App. |A.l 
read in the high-energy limit (k,j > 0): 

B...(p,m Q ,mi) ~ B...(p,0,0), (26) 
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Coo-(Pi> (Pi -Pz) ,P2,m l} m 0} m2) ~ Coo-(0, (pi -p 2 ) ,0,0,0,0), (27) 

(_l)*+i(A; - - 1)1 

C 1 ...i2...2(pi,P2,m ,m h m 2 ) -, ■ ,' , ^ ' , (28) 

(pi -p 2 ) 2 (k + j)\ 

k ] 

O x . . .i(pi,p 2 ,m ,mi,m 2 ) ~ t — --{ J Bi...i(pi -p 2 ,0,0) - £i...i(pi,rao,rai)}. (29) 

V (P1-P2) 2 



The most important functions are explicitly given by: 



and 



S (P, 0, 0) = A - In (-^T^) + 2 ' ( 30 ) 

£i(p,0,0) = _ls (p,0,0) } (31) 

Sn(p,0,0) = ls (p,0,0) + l (32) 

2 2 

B 00 (p, 0, 0) = -^S (p, 0, 0) - ^, (33) 



Ci2{pi,P2,mo,m 1 ,m 2 ) ~ — r=, (34) 

2 (Pi -P2) 2 

Cii2(pi,P2,wo,mi,rn 2 ) ~ -— (35) 

6(pi -P2) 2 

Coo(Pi,P2,rn ,m 1 ,m 2 ) ~ -5 (pi - p 2 , 0, 0) + -, (36) 

1 1 

C 001 (p 1 ,p2,m ,m 1 ,m2) — £? (pi -p 2 ,0,0) - — . (37) 



3 High-energy approximations of one-loop 4-point functions 



The definitions of the scalar and tensor 4-point functions are given in App. [A.l| . The 



high-energy expansion of the scalar 4-point function can also be evaluated using the 
Mellin-transform technique. Therefore, we sketch only the most important aspects and, 
in particular, the differences to the case of the 3-point function discussed above. For the 
tensor 4-point functions the Mellin-transform technique breaks down and we have to rely 
on tensor-integral reduction. 

3.1 Calculation 

High-energy approximation of the scalar 4-point function means that both \{k\ + A; 2 ) 2 | 
and \{k\ + ^4) 2 1 are large compared with the internal and external masses: 

I (Aii + fc 2 ) 2 |, \(k x + k A ) 2 \ > m 2 , ml, ml, m 2 3 , k\,k\,k\,k\. (38) 

Therefore, we need a double Mellin transform in order to evaluate the expansion. The 
scalar 4-point function has the Feynman-parameter representation 

D (ki, ki + k 2 , h + k 2 + k 3 , m , m x ,m 2 , m 3 ) = 

f°° A A A A 8(1-X -Xi-X2-X 3 ) 

/ dx dxidx2dx 3r - r — — — 39 

Jo {g{x , Xi, x 2 , x 3 ) - (0-1 + ie)rix x 2 - (cr 2 + ie)r 2 XiX 3 y 



with 

3 

g(x , X!,x 2 , x 3 ) = ^2 m n x n - k 2 x xi - k\x\x 2 - k\x 2 x 3 - k\x 3 x G - ie, (40) 

n=0 

n = \{ki + k 2 ) 2 \, r 2 = \(ki + k A ) 2 \, ox = (h + k 2 ) 2 /ri, a 2 = (k x + k 4 ) 2 /r 2 , where 
fc 4 = — ki — k 2 — k 3 . As above, we have included an infinitesimal imaginary part in 
the terms involving the large parameters (k x + k 2 ) 2 and {ki + /c 4 ) 2 . Thus, we can use the 
generalized Feynman-parameter representation ( |10"D in order to perform the double Mellin 
transform resulting in: 



A<[A>,£i,6 



r(6)r(6)r(2-6-6) 



(-cri-ief (-<x 2 -ie) 6 

f°° , , , , 5(1 -x -xi -x 2 -x 3 ) 
x / dxodxidx2dx3 7 — 7 — 7 — t-- (41) 

"'O ^(^0) ^1) x 2i X 3 ) 2 ~^~&XqXi X 2 X 3 



From the explicit results of Ref. [ |2| we expect the following form for the high-energy 
expansion of the scalar 4-point function: 



2 a 



Do ^^(riM^ + O - +0 - . ri,r 2 -oo (42) 

nf2 ^oJij! V r i r 2/ \ r i r 2/ 

with C22 = C12 = C21 = 0. Accordingly, the Mellin transform should have the following 
meromorphic structure at £1 ~ 1 and £2 ~ 1- 



M [D , a, 6] = £ (1 _ _ + - £0°) + - 6)°). (43) 



Note that terms that involve only poles in one of the variables are irrelevant for the 
asymptotic expansion, as they do not contribute to the inverse Mellin transform. 

In contrast to the scalar 3-point function, we divide the integration domain of the 
Mellin transform into four sectors and expand each sector separately: 

X > X U X 2 , X 3 , Xx > x 2 , x 3 , Xq, x 2 > x 3 , x , x u x 3 > x , x u x 2 . (44) 

The contribution of each sector can be associated to one side of the box and must be 
decomposed into several cases like for the scalar 3-point function. The third Gamma 
function in fl4l|) involves a pole 1/(2 — £1 — £ 2 ). This pole does not fit into the structure 
flUf). It is essential for our approach that this pole is cancelled in the total result for the 
Mellin transform, although it is present in each single sector. 

Because all sectors yield integrals of the same form, we consider only the sector xq > 
xi,x 2 , x 3 . We substitute Xi — > xoXi, i = 1, 2, 3, in fl4l[) and obtain for the integral: 

So{£i,£,2) — / dxidx 2 dx 3 t — 7 — 7- (45) 



with 



h(xi,x 2 , x 3 ) = ml + (m 2 + m 2 Q - k\)xi + (m\ + m 2 Q )x 2 + (m 2 + m 2 Q - kj)x 3 

+ (m\xi + m\x 2 + m\x 3 )(x\ + x 2 + x 3 ) - k\x\x 2 - k\x 2 x 3 - ie. (46) 
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As for the scalar 3-point function, we consider first the case wig 7^ and use 

1 111 

+ 



h(x 1 ,x 2 ,x 3 ) 2 -^-& h(x u 0, 0) 2 ~£i-& h(0,0,x 3 ) 2 -^-& h(0,0,0) 2 -^-b 

+ (2 - ^ - 6)(0(x 2 ) + 0(six 3 )). (47) 

The terms that have not been written down explicitly involve only poles in one of the 
variables 6 and are hence not relevant for the asymptotic expansion. Thus, we find 



;2-6-6) f 1 



1 -eoa- 6) U(o,o, 0)1-6(1 -e 2 ) 2 



x v Mo, o,o)y 1/0 x V Mo, 0,0) y J 
1 (2 - 6 - - 6)°) + - 6)°)}- (4* 



1 3 



The factor (2— 6— 6) in the first and third terms cancels the pole contained in r(2— £1— £2) 
in (£11]). If we combine all four sectors this pole is also cancelled in the contributions of 
the second term: 
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2(2-6-6) f 1 + 1 1 



(1-6X1-6) l(i-6) 2 (i-6) 2 (i-6)(i-e 2 )j 

+ 0(2-6-6). (49) 

This holds as well in the other cases with m 2 , = considered below. 

It is then straightforward to construct the expansion of the Mellin transform and eval- 
uate the high-energy approximation by comparing fl42 ) and (£E|). Because the calculation 
is lengthy but analogous to the one for the scalar 3-point function, we only make some 
remarks. The integrals involving the logarithms in fl48f) lead to contributions of the form 



I D (kl m , rm) = - I' dx- In ( 1 + <±<A X + -^— x 2 \ . (50) 
Jo x \ rriQ — ie ttiq — ie J 

Each of these integrals can be associated with one corner of the box and for each corner 
two integrals appear in the sum of all sectors. Using the relation 

7T 2 

I D (kl, m , mi) + I D (kl, m u m ) = I c (k 2 , m , mi) + I c (k\, m x , m ) - y (51) 

with Wq, m 2 7^ 0, all these integrals can be expressed by those appearing in the asymptotic 
expansion of the 3-point functions. Finally the asymptotic expansion of the scalar 4-point 
function can be written in terms of the asymptotic expansions of the four corresponding 
scalar 3-point functions and a universal term. 

Next we consider the case m 2 , = 0. We decompose the integration domain of the 
integral S into three sectors 

x 1 >x 2} x 3 , x 2 >x 3 ,x 1 , x 3 >x 1 ,x 2 , (52) 



9 



and arrive at 



So = dxido^dxa 
Jo 



h\{x\, x 2 , x 3 ) 2 ~^~^ 2 x & x ^ a;| 



■i 



with the polynomials 



f 1 1 
+ / dxidx 2 dx 3 7 — t — f 

Jo h 2 (xi,x 2 ,x 3 ) 2 ^^^i 2 x\ 2 x 2 x 3 
+ / dxidx2d^3 t — t — 7- (53) 

Jo h 3 (xi,X 2 ,X 3 ) 2 ~^~^ 2 x\ 2 X 2 X 3 



h 1 (xi,x 2 ,x 3 ) = h(xi,x 1 x 2 ,x 1 x 3 )/x 1 , 
h 2 (xi,x 2 ,x 3 ) = h(x 2 x 1 ,x 2} x 2 x 3 )/x 2 , 

h 3 (x 1 ,x 2 ,x 3 ) = h(x 3 x 1 ,x 3 x 2 ,x 3 )/x 3 . (54) 



If m\ k\ and m 3 ^ k 2 , the first and third integral are analogous to (45). The second 
integral can be evaluated by performing an expansion of the integrand with respect to 

(2-6-6), 

. 1 da;idx 2 da;3 > f > = 7 ^-tt; + (2 - 6 - f 2 )0((l - 6)°), (55) 

^2(^1,^2, a;3) 2- «~«a;f x| 2 x| 2 (1 - 62) 

which is true no matter whether h 2 ( 0,0,0) = or not. Notice that (f49|) holds also after 
performing sector decomposition, because the pole 1/(1 — £ 2 ) 3 °f the second integral in 
(|53|) [see (|55|)1 cancels one of the poles in the first and third integral [see ([Hp ]. 



If m\ = k\ 7^ or m| = 7^ we obviously have to apply a further sector decompo- 
sition to the first or third integral, respectively. In this way, all 4-point functions which 
contain neither IR nor mass singularities can be approximated for high energies. Again 
the singular ones are covered by our approach once they are regularized with masses. It 
turns out, that in all cases the scalar 4-point function can be expressed in a universal way 
in terms of the four corresponding 3-point functions. 

In contrast to the case of the scalar 4-point function, the pole 1/(2 — 6 — £ 2 ) does 
not cancel in the Mellin transform of the coefficient functions of the tensor 4-point inte- 
grals. Because the poles of the variable £1 and £ 2 are coupled in the Mellin transform, 
the resulting asymptotic expansion depends on the order of the inversion of the Mellin 
transforms with respect to £1 and £ 2 - Integrating first over £1 and then over £ 2 results 
in an approximation for f] > r 2 > m 2 ,k 2 , integrating first over £ 2 and then over 6 
yields an approximation for r 2 ^> r\ ^> m 2 ,k 2 . While these approximations coincide 
for the scalar 4-point function, where the poles in £1 and £ 2 factorize, they differ for the 
tensor 4-point functions. Therefore, we have evaluated the approximations for the tensor 
4-point integrals by reducing them to scalar integrals with FeynCalc [ 0] and determining 
the high-energy limit by using the approximations for the scalar integrals with the help of 
Mathematica [12]. As can be seen from the results in Sect. 3.2.2| , the pole 1/(2 — 6 — £ 2 ) 



in the Mellin transform corresponds to a factor \/{a\r\ + cr 2 r 2 ) (1/u in the notation of 



Sect. p.2.2| ) in the high-energy approximation of the integrals. 
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3.2 Results 

The results are valid for 

\s\, \t\, \u\ S> k\, k\, k\, k\, m\, m\, m\, m\ (56) 

with s — (ki + ki) 2 , t — (ki + k 2 ) 2 and u — (ki + k 3 ) 2 . All 3-point functions in this section 
are implicitly understood in the high-energy limit. 

3.2.1 Scalar 4-point function 

The high-energy approximation of the scalar 4-point function can be reduced to the 
approximations of the four corresponding 3-point functions as follows: 



D (h, ki + k 2 , h + k 2 + k 3 , mo, mi, m 2 , m 3 ) ~ — - 

st 



—s — ie 



+ ^Coih, -k 2 , mi, mo, m 2 ) + jC (k 2 , -k 3 , m 2 , m x , m 3 ) 

+ -C (k 3 , -h, m 3 , m 2 , m ) + yC (/c 4 , ~h, m , m 3 , mi). (57) 

Note that this result holds for all non-singular cases, in particular for arbitrary masses. 
The necessary split-up into different cases has to be done only at the level of 3-point 
functions. Specific examples of this general result can be found in App. A of the second 
paper of Ref. [ |J. 

3.2.2 Tensor 4-point functions 

The high-energy approximations of the coefficient functions of the tensor 4-point in- 
tegrals, defined in App. [A.l|, can be written as: 



D 1 ^~(L 3 +^)--C 0i (58) 
2su s 

D n ~ -t^(L 2 + it 2 ) — —L + -(C + C x + C 2 ), (59) 
2su A su s 



D 12 ~ —-^—(L 2 + vr 2 ) + —L- -C 2 , (60) 
2u A su 

1 ,,n . 1 



D 13 ~ _(L 2 + 7r 2 )-— L, (61) 
2u z su 

t 2 ,t2 2^ 3t + S r t + 2s 

Dm o(L + it) ttL H 

2sm sV ; 2su 2 2tsu 

- -(Co + Ci + 01-^11-^22), (62) 
s 

Cll2 .-^ + ^) + ^£L-J--i C22 , (63) 



Dm ~ ^(I 2 + * 2 ) - \-^L - -L, (64) 
2ir 2su z 2su 



^123 ~ ^(L 2 + 7T 2 ) _ -Li + J_, (65) 
4W 3 M z 2sm 
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t 3 , r2 2N lit 2 + Its + 2s 2 r 2t 2 + 7ts + 4s 2 
-Dim ~ i(L + 7T ) 5 D H 

+ -{C + C l + C 2 - C n - C 22 + C in + C 222 ), (66) 
s 

n (T 2 , 2x , 2t 2 -5ts-s 2 r 2t + S 1 

£>iii2 — (L+vr) + — — L-— — C 222 , (67) 

2m 4 bsu 6 2tu z s 

Ou U ~ S? (£» + «') 6i^ L "2^' (68) 

ts /r9 9x t 2 + 5ts-2s 2 t 2 -ts + s 2 

^1122 ~ j(L + 7r) H L 

2u 4 osir 6tsu 2 

+ -(C 22 + C 222 ), (69) 

^1123 — ^ — (d + 7T 2 ) - ~rir + ^r^' ( 7 °) 

6m 6m j 6sm 

^nsa ~ + - 2 ) - - ^ (71) 

6m 4 osir 6su z 

s 2 -2ts. r2 2 , £-5s r t + 4s 

-D1223 ~ 3 — (L + 7T ) H — L H — . (72) 

1223 6m 4 V ; 6m 3 6sw 2 ' V ; 

D 00 ~--^(L 2 + 7r 2 ), (73) 
4u 

D mi ~-J-{L 2 + ^)-^-L, (74) 

4m 

fl0011 „ + A -^L-p ,75) 

°™~^»+*)-h L -h (77) 

Cmm „ __|_ (L2 + ^)-^ + | + 1*«,0,0), (78) 

where 

The arguments of the 3- and 4-point functions are given by: 

D... = D...(h, h + k 2 , ki + k 2 + k 3 ,m ,m 1 ,m 2 ,m 3 ), 

C... = C...(ki,—k 2 ,mi,mo,m 2 ). (80) 

For the 3-point integrals the approximations given in Sect. 1273 have to be inserted. Note 



that only one kind of 3-point functions appears. Moreover, the coefficient functions with 
indices 00 or 1 and 3 are independent of all internal and external masses and do not 
involve dilogarithms. 

4 Conclusion 

We have calculated high-energy approximations of the basic one-loop integrals necess- 
ary for scattering processes of two on-shell particles into two on-shell particles (2 — > 2 
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processes). All kinematical variables have been assumed to be large compared with the 
internal and external masses, and all terms with the leading power of the energy have 
been kept. 

We used the Mellin-transform technique and overcame the difficulties in expanding the 
Mellin transform by sector decomposition. As a consequence, the results for the scalar 
integrals fall into several different cases. It turned out that the high-energy approximation 
of the scalar 4-point functions is just given by the four corresponding 3-point functions 
and an universal term. The tensor 2- and 3-point integrals were expanded with the 
help of the Mellin-transform technique and checked by using tensor-integral reduction. 
The approximations for these integrals are, in contrast to the one for the scalar 3-point 
function, very short and include no dilogarithms. The tensor 4-point integrals were first 
reduced to scalar functions and then approximated resulting in compact expressions in 
terms of 3-point functions and logarithms. 

We have listed high-energy approximations for the complete set of scalar one-loop 
integrals that appears in 2 — > 2 processes. In addition, we have provided explicit results for 
the corresponding tensor integrals. In the case of the 4-point function we have restricted 
ourselves to tensor integrals with at most four Lorentz indices. This set is sufficient for 
calculations in the 't Hooft-Feynman gauge. The tensor integrals with more indices that 
appear in more general gauges can be easily reduced to those given here using tensor- 
integral reduction. 

With the help of our results and tensor-integral reduction any one-loop matrix element 
for 2^2 scattering processes can be approximated for energies much larger than the 
internal and external masses. Moreover, our approximations are applicable to the decay 
of a very heavy particle into two light ones when all internal masses are also light. Since 
our results contain only the leading terms they allow to extract the leading contributions 
to matrix elements in the high-energy limit, i.e. the non- vanishing terms for 2 — > 2 matrix 
elements. 



A Appendix 

A.l Definition of one-loop integrals and tensor-coefficient functions 

We use dimensional regularization denoting the space-time dimension by D. The 
one-loop functions are defined by: 

Aoh) = / d D q ±- (Al) 

m z J D 

B { o,^u}(pi,mo, mi ) = — J d q — rg-g- , (A2) 

C { o,n,^,nup}{pi:P2,mo,m 1 ,m 2 ) = — / d q , (A3) 

17T J L>qL>iL>2 

,p,pu, pup, pupa} 

(pi,P2,P3,mo,m 1 ,m 2 ,m 3 ) = 

(2vr / u) 4 - D /-^ {l, %, q^, q^qp, g^gpgq} 



ivr 2 J DqD^Ds 
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with D = q 2 — TOq + ie, Di = (q + pi) 2 — to| + ie, i > 1. The tensor integrals can be 
decomposed into Lorentz tensors constructed of (linearly independent) external momenta 
Pip and the metric tensor g^ u and tensor-coefficient functions as follows: 

Bp = B lPlfl , (A5) 
By. v = B u p lfl p lu + B 00 g^, (A6) 

C, = j2C lPl „ (A7) 
i=i 

2 

Cfiu = CijPif+Pjv + Coog^, (A8) 

2 2 

^Wp = 51 C i j k pi fJ pj V p k p + ^CQQi{g tlv pi P + gu P Pifj, + g^pPw}, (A9) 

i,j,k=l i=l 
3 

D M = ^ APv, (A10) 
i=i 

3 

= H Cijp^pjy + D^g^, (All) 

3 3 

^Vp = 51 DijkVi^VjuVkp + D Q0i{9puPip + 9v P Vip + 9ii P Viv), ( A12 ) 

i,j,k=l i=l 
3 

PkpPia + D QQ0Q {g^ u g pa + + 9^9up} 

i,j,k,l=l 
3 

+ 5Z D mj {g^p ip p ja + gppPi U p ja + gpaPiuPjp 

+ 9upPipPja + guaPipPjp + 9pcrPinPjv}- ( Al3 ) 

Since the tensor-coefficient functions are scalars and depend only on kinematical in- 
variants, internal and external masses, we can write their arguments as: 

C. ..(p 1 ,p 2 ,m ,m 1 ,m 2 ) = C...{p\, (p 2 - Pi) 2 ,pl, mi, m , to 2 ), (A14) 

D. ..(ki,h + fc 2 , &i + /c 2 + fc 3 , mo, mi, m 2 , m 3 ) = 

= D...(k 2 ,kl,k%,kl, (ki + k 2 ) 2 , (fa + k 4 ) 2 ,m ,m 1 ,m 2 ,m 3 ), (A15) 

where for convenience we introduce the momenta of the external particles (cf. Fig. |I]). 

A. 2 Feynman-parameter representation of one-loop functions 

The UV divergences of the following integrals are contained in the quantity 

A=£-^- 7j5 + lii(47r), (A16) 

where 7^ is Euler's constant. 
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q + ki + k 2 + k. 



Figure 1: Conventions for 4-point integrals 



A. 2.1 1- and 2-point functions 

The tensor-coefficient functions of the 1- and 2-point integrals can be represented as 
(i > 0): 



A (m )=ml A-ln P^l+l 



Bi... 1 (pi,m ,m 1 ) 



i + 1 



A — / dx dxix\5(l — x — X\) In 



(A17) 



, (A18) 



■ iii TflrX 

Sooi..^(pi,mo,mi) = (-1)*^ | — , n\ + C , \\ ~ C , n\C \ q\ 

s -v- / 2ll(j + l)(i + 2) (* + 2) + 2)0 + 3) 



(A + l) 



dxodxia; i 1 5(l — xq — X\)M% In 



(A19) 



with = m^xo + m^i — p^xoXi — ie. Note that the scalar 2-point function Bq is given 
by ( gig ) with i = 0. 



A. 2. 2 3-point functions 

The tensor-coefficient functions of the 3-point integrals C...(pi,p2,mo,mi,i7i2) allow 
for the following Feynman-parameter representations > 0): 



Ci. 



12- - 2 



-(-l)*' f dxodx.do;^ 1 - X ° " Xl - X2) 

JO 



(A20) 



Cooi--- 12- ■■ 2 — (— 



2 1 (i + j + 2)! 



d^odxida^o;^^! — xq — x\ — X2) In 



'Mc 



(A21) 
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with Mq = m^xo + m\x\ + m\x 2 — p\x§X\ — p\x$x 2 — (pi — p 2 ) 1 X\X 2 — ie. The scalar 
3-point function is given by (|A20|) with i = j = 0. 



A. 2. 3 4-point functions 

The Feynman-parameter representations for the tensor-coefficient functions of the 
4-point integrals D...(ki, k\ + k 2 , k\ + k 2 + k 3 , mo, mi, m 2 , m 3 ) read (i, j, k > 0): 



POO ,y,? ,-y,3 rpk S. ( 1 \ S *J rp \ 

,} L^3K3 { irJ '\L ^1^3 1 2 3 l M4 Ln=D nJ , (A22) 

i j k ^ 

I POO rr-i rf.k (\( 1 3 \ 

^ooi_i2_23_3 = -^(-l) 4+J+fe / o dx Q d Xl dxAx ^ W{l M2 Xn) , (A23) 



1 , f it 

MjOOOl-- 12- - 23- •• 3 — tI — i 



i j h 



4 V ' [(i + j + k + 3)\ 

r°° ( M 2 \ 1 

- ^ dxodx^dxsxlxix^il - J2l= Xn) In l-^-j j (A24) 

3 

with M|> = X) m n x n — fciXoXi — &|xi£2 — k 2 x 2 x 3 — klx 3 x — (ki+k 2 ) 2 x Q x 2 — (ki+k4) 2 xix 3 — ie 

n=0 

and /C4 = — ki — k 2 — k 3 . The scalar 4-point function is given by ( |A22|) with i = j = k = 0. 
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